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ABSTRACT 


A study has been made to determine the accu- 
racies which could be expected from onboard primary 
navigation systems for the lunar module following 
the concentric flight plan. Emphasis was placed on 
comparing a rendezvous radar with an optical tracker 
as a source of navigational information. The com- 
parison was drawn from linearized error analyses 
of the navigation systems as they would function in 
various phases of the concentric flight plan. The 
results of the investigation indicate that there is little 
difference between rendezvous radar and optical 
tracker performance in the primary navigation system 
in terms of navigation accuracy. Furthermore, for 
either sensor, the effect of reduced accuracy and the 
effect of varying the measurement sampling period 
on navigation accuracy are of little significance. 
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A LINEARIZED ERROR ANALYSIS OF ONBOARD PRIMARY 


NAVIGATION SYSTEMS FOR THE APOLLO 


LUNAR MODULE 


By Jerrold H. Suddath, Robert H. Kidd in, 
and Arnold G. Reinhold 
Manned Spacecraft Center 


SUMMARY 


A study has been made to determine the accuracies which could be ejected from 
onboard primary navigation systems for the lunar module following the concentric flight 
plan. Emphasis was placed on comparing a rendezvous radar with an optical tracker 
as a source of navigational information. The comparison was drawn from linearized 
error analyses of the navigation systems as they would function in various phases of 
the concentric flight plan. The results of the investigation indicate that there is little 
difference between rendezvous radar and optical tracker performance in the primary 
navigation system in terms of navigation accuracy. Furthermore, for either sensor, 
the effect of reduced accuracy and the effect of varying the measurement sampling 
period on navigation accuracy are of little significance. 


INTRODUCTION 


The purpose of this study was to compare a rendezvous radar with an optical 
tracker as a source of navigational information in an onboard primary navigation sys- 
tem for the Apollo lunar module (LM). The comparison was drawn from error anal- 
yses of one system with an optical tracker and another system with a rendezvous radar. 
Included in this comparison are the effects of sun and moon interference on the sen- 
sors, the inertial platform alinement accuracies of the two systems, the co mman d and 
service module (CSM) ephemeris uncertainties, and the navigation sensor reference 
misalinement. 

The error analyses were performed by using a digital computer to calculate the 
covariance matrix for the errors of a navigation system (using either the radar or 
tracker) as a function of time along the trajectories of interest. It was assumed that 
the navigation system would estimate the state (position and velocity) of the LM meas- 
ured relative to the CSM and that the estimation would be done with a Kalman filter 
(ref. 1). Results are presented in the form of time histories of root mean square (rms) 
relative position and velocity estimation uncertainties for the trajectories considered. 
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unit vector in direction of sun 
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SCOPE OF THE STUDY 


To compare a rendezvous radar with an inertially referenced, optical line- of - 
sight tracker as a source of navigational information for the LM, it was necessary to 
consider the salient characteristics of the overall LM mission. In this study, the 
LM mission was taken to be that of following the concentric flight plan (fig. 1). The 

nominal sequence of events on this flight 
plan was taken to be the following : 

(1) LM separation from CSM 

(2) Hohmann descent 

(3) Powered descent 

(4) Landing 

(5) Nominal launch 

(6) 90° burn to adjust apocynthion 

(7) Circularization 

(8) Transfer maneuver 

(9) Terminal rendezvous 

In addition to the nominal sequence 
of events, there were the following con- 
tingencies to be examined : 

(1) Abort 12 minutes after LM separation; 70° direct transfer to rendezvous 

(2) Abort 35 minutes after LM separation; 140° direct transfer to rendezvous 

(3) Abort from start of powered descent 

(4) Abort from hover 

(5) Late launch 

The necessity of having to consider all these phases of the mission arose because 
there are physical limitations on both sensors which restrict their ability to track. 

For example, the optical tracker cannot track if the angle between the lines of sight 
(LOS) to the sun and to the CSM is less than some critical value which is shown in fig- 
ure 2 as a function of relative range. Also, the tracker cannot track the CSM against a 
fully ill um inated lunar background at ranges greater than 40 nautical miles. Moreover, 
neither the tracker nor the radar can track at ranges greater than 400 nautical miles. 
Thus, every phase of the mission had to be examined to see if these limitations created 
any major problem areas. 


Sun elevation angles 

at landing site: 15° to 45* 


Circularization 

altitudes: 30 to 65 n.mi. 



Figure 1. - Concentric flight plan. 
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Angle bias Cl#) = 0.2 mr 
Random angle G# ) = 0.45 mr 


METHOD OF ANALYSIS 



0 X00 200 300 400 

Relative range, ami 


Figure 2. - Optical tracker error model 
and sun angle constraint. 


lated only in the sense of computing their 
process. 


The error analyses were performed 
by using a digital computer to calculate 
the covariance matrix for the errors of a 
navigation system (using either radar or 
tracker) as a function of time along the 
trajectories of interest. The computer 
program used is described in the appendix. 

The rms relative position uncer- 
tainty was obtained by taking the square 
root of the sum of the first three diagonal 
elements of the covariance matrix. Simi- 
larly, the rms relative velocity uncer- 
tainty was obtained by taking the square 
root of the sum of the second three diag- 
onal terms of the covariance matrix. 
Navigational measurements were simu- 
effects on the statistics of the estimation 


To do this type of error analysis, certain a priori statistics must be given. 

Also, statistical error models are necessary for the radar, tracker, accelerometers, 
coupling data unit, and inertial platform and sensor reference misalinements with 
respect to the inertial platform. The error models used in this study are discussed 
in the next section. 


ERROR MODELS 


Optical Tracker 

In using the model for the optical tracker, it was assumed that the tracker would 
measure the azimuth and elevation angles of the CSM relative to an inertial coordinate 
system (the inertial platform). The sun interference constraint is plotted in figure 2. 
Other pertinent data are 


Bias (1 cr) = 0.2 milliradian (mr) 

Noise (1 cr) = 0.15 mr 

For the tracker, specification data indicated that sensor reference misalinement 
errors were negligible. 
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Rendezvous Radar 


In using the model for the radar, it was assumed that the radar would measure 
either relative range or relative range rate and shaft and trunnion angles (defining the 
line of sight to the CSM) relative to the sensor reference. Pertinent data for the radar 
are listed below. 


Range bias (3 a) = 1500 feet 

Range noise (3 a) = greater of 0. 25 percent or 300 feet 

Range-rate bias (3 cr) = 3 ft/sec 

Range -rate noise (3 a) = greater of 0. 25 percent or 1 ft/sec 

Angle bias (3 a) = 7 mr (including sensor reference misalinement) 

Angle noise (3 cr) = 2 mr for range < 200 nautical miles linearly 

increasing to 4 mr at 400 nautical miles 

For the rendezvous radar, the sensor reference misalinement errors (specifica- 
tion values) are significant to the point where vehicle attitude must be known to repre- 
sent the errors in a measurement. This being the case, vehicle attitude profiles had 
to be assumed in the simulations of the radar. The nominal attitude profile used in 
all radar simulations is described as follows: 



(1) No roll or yaw; LM X- Z plane 
always coincident with orbital plane (fig. 3) 

(2) Pitch vehicle so as to maintain 
LM Z-axis along line of sight from LM to 
CSM 


The reasons for selecting this attitude 
profile are the following: 

(1) With the vehicle Z-axis along the 
line of sight, the radar is always bore- 
sighted down the Z-axis. This attitude 
tends to minimize the radar boresight drift 
error, since it is a function of the shaft 
and trunnion angles which, in this case, 
are nominally zero. 


Figure 3. - LM attitude profile. 
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(2) Since the normal field of view for the astronauts is in the direction of the LM 
Z-axis, this attitude profile keeps the CSM in their field of view at all times. 

(3) This attitude profile also simplifies the mechanization of the onboard naviga- 
tion system. 


Inertial Platform Errors 


When a AV maneuver is made, the errors in the estimate of the applied AV 
increase the uncertainties in the estimate of the vehicle velocity. Since the inertial 
platform measures the applied AV, the error estimating the applied AV comes from 
the inertial platform errors. Therefore, after a AV maneuver, the covariance ma- 
trix for the error in the navigational estimate must be updated according to the error 
model for the inertial platform. The data used for the inertial platform errors (spec- 
ification values) are given below. 

Accelerometer bias (1 cr) = 0. 00656 ft/sec^ 

Accelerometer scale factor (1 a) = 100 parts per million 
Gyro drift (1 cr) = 0. 5 mr/hr 
Platform alinement: Optical Tracker (1 a) = 0. 2 mr 
Platform alinement: Rendezvous Radar (1 a) = 1 mr 


The inertial platform was assumed to be realined within 15 minutes of each AV 
maneuver. This essentially eliminates the effect of gyro drift on the measurement of 
the AV maneuver compared with the effect of alinement inaccuracy. An alinement 
optical telescope is used to aline the platform in the rendezvous radar system. 


RESULTS AND DISCUSSION 


Various trajectories associated with the concentric flight plan were simulated by 
using the error models presented in the previous sections and the computer program 
described in the appendix. The three basic cases considered for each trajectory were: 

(1) Navigational information provided by optical tracker measuring azimuth and 
elevation angles 

(2) Navigational information provided by rendezvous radar measuring shaft and 
trunnion angles and range 

(3) Navigational information provided by rendezvous radar measuring shaft and 
trunnion angles and range rate 

The results of the simulations are discussed in the following sections. 
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Nominal Launch 


The geometry of the nominal launch trajectory is represented in figure 1. The 
trajectory consists of the following sequence of events : 

(1) Powered ascent from launch site (taken to be 0° long. , 0° lat. ) to ascent in- 
jection 

(2) Coast through approximately 90° central angle 

(3) At approximately 90°, a AV maneuver is made to adjust the altitude and 
longitude of apocynthion 

(4) Circularization of the orbit at apocynthion and coast to transfer 

(5) At the transfer point, a AV maneuver is made to put the LM on a 140° 
direct transfer to rendezvous 

The simulation of this phase was initiated at ascent injection. Navigational meas- 
urements were assumed to be made every 3 minutes if the constraints on the sensors 
were satisfied. 

Results of the simulations are shown in figure 4. Figure 4(a) is a semilog plot 
of the rms relative position uncertainties in feet versus time after ascent injection in 
minutes. 

The curve at the top of the figure is for the case where no measurements are 
made throughout the trajectory. In this case, the rms relative position uncertainty 
begins at about 3700 feet and simply grows to about 68 000 feet at the time of rendez- 
vous. 


The curve at the bottom of the figure is for the case where the optical tracker 
was used to make measurements. The initial uncertainty for this case is about 400 feet 
less than for the other cases plotted on the figure, which results because prior to 
launch the platform is realined, and alinement with the optical tracker system is more 
accurate than alinement with the radar system. Therefore, when the alinement errors 
are propagated through the powered ascent, the tracker system has smaller injection 
errors than the radar system. 

The optical tracker experiences sun interference, which is indicated in figure 4(a) 
as occurring at 67 minutes and lasting until 78 minutes after injection. This trajectory 
was computed assuming that the sun elevation angle (relative to the lunar local hori- 
zontal) at the landing point was 45°, and the lunar stay time was 36 hours. However, 
other trajectories were computed with the elevation angle being varied upward from 
an elevation of 15°, and no significant difference was noted in the performance of the 
tracker system. Therefore, this case was presented as being typical of any sun eleva- 
tion angle currently being considered for the concentric flight plan, and it was con- 
cluded that sun interference with the optical tracker is not a problem on the nominal 
launch. 
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Rms velocity error, ft/sec Rms position error, ft 


100 000 


The two curves in the middle of fig- 
ure 4(a) are for the two cases where the ren- 
dezvous radar was the sensor. In both cases, 
shaft and trunnion angles were measured so 
that the cases differ only in that range was the 
third measurable in one case, whereas range 
rate was the third measurable in the other 
case. The two types of radar provided very 
little difference in the navigational accuracy. 

Figure 4(b) is a semilog plot of the rms 
relative velocity uncertainties in feet per sec- 
ond versus time after injection in minutes. 

The description of this plot is virtually the 
same as the preceding description of fig- 
ure 4(a). 

The "lumps" on the radar curves are 
caused by out -of -plane errors and are easily 
seen by examination of figure 4(c) which is a 
plot of the components of the relative position 
uncertainties along the line of sight, out of 
plane, and in the orbital plane normal to the 





(c) Trajectory error with no measurements- nontin a l launch. 


Figure 4. - Relative position and velocity error for nominal launch: 
range radar (RR), range rate radar (ftR), optical tracker (OT). 
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line of sight. Although the plot is for the case where no measurements are made, the 
periodic character of the out -of -plane error is very evident in the radar errors plotted 
in figure 4(a). 

In the navigation system of this study, the relatively large radar out -of -plane er- 
rors were due to a misalinement bias. If necessary, the error could be reduced sub- 
stantially by simply estimating this bias, which would essentially eliminate any 
difference between the radar and tracker accuracies, but would impose some additional 
computational requirements on the radar system. Overall comparison of the curves on 
figures 4(a) and 4(b) shows that there is little difference in navigation accuracy using 
either the optical tracker or the rendezvous radar in the primary navigation mode. 


Late Launch 

The geometry and sequence of events 
for the late launch trajectory are basically 
the same as for the nominal launch trajec- 
tory. The only substantial difference is the 
position of the CSM at launch. Figure 5 is 
a semilog plot of the rms relative position 
uncertainty in feet versus time after injec- 
tion in minutes. The description of this 
plot is similar to the description of fig- 
ure 4(a). 

During the first 21 minutes of the 
late launch trajectory, neither the optical 
tracker nor the radar can make measure- 
ments, because the relative range diming 
this period is greater than 400 nautical 
miles. Once again the optical tracker ex- 
perienced sun interference, but it was 
found that interference is not a problem 
on the late launch. 



Figure 5. - Relative position error for 
late launch: .range radar (RR), range 
rate radar (RR), optical tracker (OT). 
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Abort Trajectory Results 


The geometry of the abort from start of powered descent is shown in figure 6. A 
semilog plot of the relative position error for an abort 12 minutes after LM-CSM sepa- 
ration, which is a direct transfer through 70°, is shown in figure 7. 



Figure 6. - LM-CSM relative inertial 
position for abort from start of 
powered descent. 



Abort time from separation, min 


Figure 7. - Relative position error for an 
abort 12 minutes after LM-CSM sepa- 
ration: range radar (RR), range rate 

radar (RR), optical tracker (OT). 


The relative position error for an abort 35 minutes after LM-CSM separation is shown 
in figure 8 using a direct transfer through 140°. This plot is virtually the same as the 
other abort, with no sun interference on the tracker for this and the above abort case 
for a sun elevation at the landing site corresponding to 15° to 45°. 

Figure 9 shows relative position error for an abort from the start of powered 
descent, with the only difference about this plot being the moon interference toward the 
end of the trajectory. The error in the optical tracker increases slowly during the pe- 
riod of moon interference, but once track is reestablished, the error comes back down. 
Therefore, moon interference with the tracker does not appear to be a problem. 

The relatively large buildup of errors during the first portion of the trajectories 
is caused by the optical tracker being hampered by the existence of a low line-of-sight 
rate during the first part of the trajectory. Whenever this is the case, the tracker 


11 





Figure 8. - Relative position error for an 
abort 35 minutes after LM-CSM sepa- 
ration: range radar (RR), range rate 

radar (RR), optical tracker (OT). 



Time from separation, min 

Figure 9. - Relative position error for 
abort from start of powered descent 
range radar (RR), range rate radar 
(RR), optical tracker (OT). 


does not provide much information, and the error tends to grow. The radar error 
grows at the beginning of the trajectory because of the rather large initial bias in the 
range. Even though the bias is being estimated, a certain amount of information com- 
ing in is required to estimate this bias. Therefore, during the first part of the trajec- 
tory, most of the information is going into the estimation of bias and very little into the 
estimation of the state. Once the bias estimate becomes fairly good, the radar errors 
tend to decrease until the point is reached in the trajectory where relative range be- 
comes quite large. Since the noise in the radar range measurement is a function of the 
range, the noise is increasing and, therefore, causing errors to increase during the 
middle portion of the trajectory. 

* Little difference was noted between the navigation accuracies obtained with radar 
and optical tracker for the abort trajectories considered, which is consistent with the 
statement made earlier for the nominal launch trajectory. 
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Effect of the Sampling Period 

The effect of the sampling period on relative position error for the optical tracker 
for the nominal launch trajectory is shown in figure 10. The solid curve shows the er- 
ror if one measurement is made every 3 minutes, and the dashed curve shows the error 
if there is only one measurement made every 9 minutes. After a measurement is made 
on the 9-minute, sampling-period case, the error is very close to the 3-minute 
sampling-period curve. When this fact is remembered and the estimate only noted im- 
mediately after making a measurement, then increasing the sample period without sig- 
nificantly degrading the accuracy of the system can be concluded. 

The effect of the sampling period on relative position error for the range radar 
shown in figure 11 is very similar to the effect of the sampling period on the optical 
tracker (compare fig. 11 with fig. 10), and the same conclusions apply. 



Time from injection, min 


Figure 10. - Effect of sampling period on 
optical tracker relative position error 
for nominal launch. Sampling period 
equals time between measurements. 



Time from injection, min 


Figure 11. - Effect of sampling period on 
range radar relative position error 
for nominal launch. Sampling period 
equals time between measurements. 
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Effect of Sensor Accuracy 

The effect of optical tracker accuracy on the relative position errors on the nom- 
inal launch trajectory is shown in figure 12. The dashed curve is the performance 
obtained if the tracker has the specification value for the noise. The solid curve is the 
error if the errors were three times the specification values. While there is some 
degradation in the performance of the system with the errors being increased by the 
factor of three, the degradation does not appear significant. 

The effect of range radar accuracy on relative position errors for nominal launch 
is shown in figure 13. The dashed curve is for the radar with the specification values. 
The solid curve is for three times the specification values. The dotted curve is for 
specification noise, but the angle bias is reduced to 1 milliradian. Once again, some 
degradation occurs when the noise and bias are three times the specification values; 
but since the errors are probably significant only at the end of the trajectory, the deg- 
radation does not seem to be very serious. 



Time from injection, min 


Figure 12. - Effect of optical tracker ac- 
curacy on relative position error for 
nominal launch. 



Figure 13. - Effect of range radar ac- 
curacy on relative position error for 
nominal launch. 
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CONCLUDING REMARKS 


A study was made to determine the accuracies which could be expected from on- 
board primary navigation systems for the lunar module following the concentric flight 
plan. Emphasis was placed on comparing a rendezvous radar with an optical tracker as 
a source of navigatio nal information. The comparison was drawn from error analyses 
of the navigation systems as they functioned in various phases of the concentric flight 
plan. The results of the investigation indicate little difference in terms of navigation 
accuracy between rendezvous radar and optical tracker performance for specification 
error models in the primary navigation system. Sun interference has little effect on 
the accuracy of the primary system with the optical tracker. Furthermore, the effects 
of measurement sampling period and degraded systems have little significance. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, January 20, 1967 
914-30-50-01-72 
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APPENDIX 


The purpose of this appendix is to describe a computer program which provides 
the capability to perform error analyses of LM-CSM rendezvous trajectories. Partic- 
ular attention has been given to the simulation of the concentric flight plan, but the 
basic program has the capability of simulating virtually any type of rendezvous tra- 
jectory with only moderate alterations. 

Fof convenience, this appendix is divided into two main sections. The first sec- 
tion is a verbal description of what the program does in terms of what is computed and 
hpw it is computed. The second section presents the appropriate equations which cor- 
respond to the description given in the first section. 


Description of the Program 

Dynamics . - The CSM is assumed to be in an orbit about the moon. The equations 
of motion for the CSM are written relative to an inertial coordinate system centered in 
the moon, which is assumed to have an inverse square, central force field. This 
representation of the potential of the moon is adequate for this type of investigation. 

The LM is assumed to be on a trajectory which would lead to rendezvous with the 
CSM. The motion of the LM is described by its position and velocity measured relative 
to that of the CSM. 

An LM guidance maneuver is represented by simply adding a AV to the LM 
velocity vector. 

Navigation system . - The navigation system is simulated by assuming that meas- 
urements would be processed by a Kalman filter. The use of this filter entails lineari- 
zation about nominal trajectories, white Gaussian noise representation of random input 
variables, and a certain amount of given, a priori statistics. 

The program will compute the covariance matrix of the error in the estimate of 
the deviations from the nominal relative state of the LM. The corresponding covari- 
ance matrix for the CSM is also computed. 

Measurements . - Either optical sightings or radar measurements can be simu- 
lated. hi the case of optical sightings, a measurement consists of the simultaneous 
reading of two angles plus bias and random errors. Radar measurements consist of the 
simultaneous reading of two angles, range or range rate, plus bias and random errors. 

The angle measurements programed for the optical tracker are azimuth and ele- 
vation, as opposed to two gimbal angles for the radar shaft and trunnion. In addition, 
the angular bias errors are assumed to be initially independent, and the noise errors 
perpetually independent. This situation is only approximately correct because of the 
actual angles measured, the inertial platform misalinement, and the sensor reference 
misalinement. The drift of the gyros is programed to affect the two biases independ- 
ently, and drifts between integration steps are assumed to be independent. Additionally, 
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before a AV maneuver, the inertial platform is realined, and the bias estimation is 
reinitialized. 

Measurements are simulated only in the sense of computing their effects on the 
statistics of the estimation process. 

Me asurement constraint . - In simulating measurements made by an optical 
tracker, it should be noted that the tracker cannot track if the angle between the lines 
of sight to the sun and to the vehicle is less than some critical value. Moreoever, this 
critical angle is a function of the range between the vehicles. This effect is simulated 
in the program by prohibiting optical measurements when a constraint is not satisfied. 


Equations for the Program 

Dynamic equations. - The equations of motion for the CSM are given by 



( 1 ) 


where 


R^, = radius vector of CSM 

r c =ll s cll 

jn = gravitational constant for moon 


Let R t denote the radius vector of LM and define the relative vector 


p=r l 



( 2 ) 


The differential equation for p is 



( 3 ) 
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where 


r L~ R l||“|| p + R C 


Nominal trajectories (solutions of equations (1) and (3) for R^(t)- and p( t)) are needed 
to compute coefficients in the matrix differential equations given in the next section. 

Transition matrices . - The transition matrix for CSM deviations is denoted by 4 > 
and computed from 


where 


and 




F(t) = 


0 

a*(t) 


I 

0 


a *(t) = w 


r C^ 


— §-R c (t)R c T (t) 

r C (t) 


(4) 


where I denotes the identity matrix. The concepts of transition and fundamental 
matrices are treated in reference 2. The fundamental matrix for relative deviations 
is denoted by 0 and computed from 


dt 0 ( t? t o) " f 0 ( t} t o)’ 0 (*0’ * 0 ) - 1 

where 

f *(t) = o 


(5) 
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and 


/3*(t) 




I ' 72 ( P + R c)( P + R c) T 

L 


Since the relative state is a function of the CSM state (see equation (3)), the additional 
transition matrix for relative deviations is denoted by r and computed from 


m r (*> ‘o) - f *<‘> ‘ 0 ) * **« * (*> ‘o) 


( 6 ) 


where 


and 


’(*<,’‘0) = ° 


g*(t) = 


|"r*(t) ol 

L 0 °J 


y*( t) = 


_ M 


r T* r L 2 ^ + RC ^ P + R ^ T ’ r C 5 R ° R ° T 


CSM Errors 

It is assumed that the estimate of the CSM state is not updated during the flight. 
Therefore, if the error in the estimate at time t^ is denoted by e^, and if is 

defined as 


p k = E ( e k i k T ) 
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then 


p M-*(Wk)V T (W*k) <’> 

where it is assumed that Pq will be given. If the estimate of the CSM state is the 
nominal trajectory, then P k is also the covariance matrix of the CSM dispersions, 
and the dispersions are the negative of the errors. 

Relative state estimation errors . - The propagation of relative state estimation 
errors is done in the following way. 

Noting that 


e k+l 


=St k + re k + T \ 


the covariance matrix of the relative state estimation error before an update is 


Mk + i = ©ir k © T + TS k T T + rp k r T + © z k r T + rz k T © T 


where the arguments of 0, r, and T are 


(jk+l* *k) 


, and 


\ - E ( ? k F k T ) 
z k - E ( f k 5 k T ) 
®k - e (\\ t ) 


u k = error in measuring a AV correction at time 
null if no correction is made 


V 


(8) 
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where 


= error in estimating relative state at time 


*k 


If 0 is partitioned thus, 


0 = 


'll 

0 21 


12 

0 22 


(9) 


then T is defined by 


T = 


e 

e 


12 
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Second, compute the gain matrix K from 


K = MB 


^BMB T + R ^" 1 


( 10 ) 


where all matrices in equation (10) are understood to have subscript k+1. The covari- 
ance matrix for measurement noise is R, and the Jacobian matrix of the measurement 
is B. Thus, if z denotes the vector valued measurement given by 


z = m(p) + v 


then 


R = 



( 11 ) 
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Finally, having computed K, the matrices n and Z are updated by 


7T = (I - KB)M 

z = (i - kb) |ez k + rp k ^<*. T 


( 12 ) 


where all matrices without subscripts are understood to have subscript k+1. The 
accuracy of the estimation process is obtained from the n matrix. 

Measurement constraint . - The angle between the line of sight to the sun and the 
CSM is determined by 


£ = cos 



where u gL is a unit vector in the direction of the sun. The critical angle £ * is 


C* = 5° + 25° j^l - e “' O154 0 40 )J 

where 1 1 p 1 1 is given in nautical miles. 

The criteria which must be met for an optical sighting to be made are as follows: 

(1) K || p| | <40 nautical miles, make a measurement only if £ > 5° . 

(2) If Up’ll > 40 nautical miles, make a measurement only if £ > £*. 



REFERENCES 


1. Kalman, R. E. ; Englar, T. S. ; and Bucy, R. S. : Fundamental Study of Adaptive 

Control Systems. ASD-TR-61-27, Vol. I. Flight Control Laboratory, ASD, 
Air Force Systems Command, 1962. 

2. Struble, Raimond A. : Nonlinear Differential Equations. McGraw-Hill Book Co. , 

Inc. , 1962. 


NASA-Langley, 1967 21 S-134 


23 



"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
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